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SUMMARY 


This paper presents stress-intensity factors calculated with a three- 
dimensionaly finite-element analysis for shallow and deep semielliptical sur- 
face cracks in finite elastic isotropic plates subjected to tension or bending 
loads- A wide range of configuration parameters was investigated. The ratio 
of crack depth to plate thickness ranged from 0.2 to 0.8 and the ratio of crack 
depth to crack length ranged from 0.2 to 2.0. The effects of plate width on 
stress-intensity variations along the crack front were also investigated. 

A wide-range equation for stress-intensity factors along the crack front 
as a function of crack depths crack length, plate thickness, and plate width 
was developed for tension and bending loads. The equation was used to predict 
patterns of surface-crack growth under tension or bending fatigue loads. A mod- 
ified form of the equation was also used to correlate surface-crack fracture 
data for a brittle epoxy material within ±10 percent for a wide range of crack 
shapes and crack sizes. 


INTRODUCTION 

Surface cracks are common flaws in many structural components. Accurate 
stress analyses of these surface-cracked components are needed for reliable pre- 
diction of their crack-growth rates and fracture strengths. However, because 
of the complexities of such problems, exact solutions are not available. Inves- 
tigators have used experimental or approximate analytical methods to obtain 
stress-intensity factors for surface cracks under tension or bending loads. For 
a semielliptical surface crack in a plate of finite thickness (fig. 1), Smith 
and Alavi (ref. 1 ) , Smith and Sorensen (ref. 2) , and Kobayashi et al. (ref. 3) 
used the alternating method to obtain the stress-intensity factor variations 
along the crack front for various crack shapes. Kathiresan (ref. 4) used the 
finite-element method to obtain the same information- However, there is consid- 
erable disagreement in the values from these solutions. 

This , paper presents stress-intensity factors calculated from a three- 
dimensional, finite-element analysis for shallow and deep semielliptical surface 
cracks in finite elastic plates subjected to tension or bending loads. Details 
of the finite-element analysis are discussed in references 5 and 6, in which 
limited stress-intensity factor results for tension loads are also presented. 

The present study covers a wide range of configuration parameters. The ratio 
of crack depth to plate thickness ranged from 0.2 to 0.8; the ratio of crack 
depth to crack length ranged from 0.2 to 2.0. The effect of plate width on 
stress-intensity factor variations along the crack front was also investigated. 
The stress-intensity factors were calculated by using a nodal-force method 
(refs. 6 and 7) . This method was preferred to the commonly used crack-opening- 
displacement method (ref. 8) , which requires a prior assumption of either plane 
stress or plane strain. This introduces a potential source of inaccuracy. The 
nodal-force method requires no such assumption. 



A wide-range equation for the stress-intensity factors along the crack 
front as a function of crack depth, crack length, plate thickness, and plate 
width was also developed for tension and bending loads. In appendix A, the 
equation is used to predict surf ace-crack -growth patterns under tension or bend- 
ing fatigue loads. These predicted patterns are also compared with measurements 
reported in the literature for steel, titanium alloy, and aluminum alloy mater- 
ials. In appendix B, a modified form of the stress-intensity equation is also 
used to correlate surface-crack fracture data reported in the literature for a 
brittle epoxy material. 
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SYMBOLS 

depth of surface crack, mm 
half -width of cracked plate, mm 

crack-growth coefficients (see eqs. (AT) and (A2) ) 

half-length of surface crack, mm 

stress- intensity boundary-correction factor 

half-length of cracked plate, mm 

mode I stress-intensity factor , kN/m^/^ 

elastic fracture toughness, kN/m^/^ 

applied bending moment, N-m 

elastic magnification factor (see eqs. (Bl) and (B2)) 
number of cycles 

exponent in equation for crack-growth rate 

shape factor for elliptical crack 

remote bending stress on outer fiber, 3M/bt^, Pa 

remote uniform-tension stress. Pa 

plate thickness , mm 

Cartesian coordinate system 

stress-intensity factor range, kN/m^/^ 

Poisson's ratio 

parametric angle of the ellipse, deg 
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THREE-DIMENSIONAL, FINITE-ELEMENT ANALYSIS 


A three-dimensional, finite-element analysis was used to calculate the 
mode I stress-intensity factor variations along the crack front for a surface 
crack in a finite plate. The plate was subjected to remote tension or bending 
loads and is shown in figure 1 . 

Figure 2(a) shows a typical f inite-elonent model of a plate containing a 
surface crack. Symmetry conditions on the x = 0 and the y = 0 planes allow 
attention to be restricted to one-quarter of the plate. Two types of elenents, 
isoparametric and singular, were used in combination to model the surface- 
cracked plate. The arrangement of elements on the crack plane (y = 0) is shown 
in figure 2(b) . The isoparametric elements (linear strain, denoted as I) were 
used everywhere except near the crack front. Around the crack front, eight 
singularity elements (denoted as S) in the shape of pentahedrons were arranged 
as shown in figure 2(c). The assumed displaconent distribution of the singular- 
ity elements had square-root terms, and therefore produced a singular stress 
field at the crack front. Details of the formulation of these types of elenents 
are given in references 7 and 8 and are not repeated here. The particular 
f inite-elonent models used were established with the help of a convergence study 
performed in references 5 and 6. The models had 4300 to 4800 degrees of freedom 


Loading 

Two types of loads were applied to the finite-element models of the surface 
cracked plate: remote uniform tension and remote bending. The ranote uniform- 
tension stress is in figure 3(a); the remote outer-fiber bending stress 

in figure 3(b) is calculated from the applied bending moment M. 


Stress-Intensity Factor 

Only loads which cause mode I deformations were considered. The mode I 
stress-intensity factor Kj for any point along the surface-crack front was 
taken to be 


ra a c 


Ki = Sidir - Fi. 


t c b 


( 1 ) 


where the subscript i denotes either tension load (i = t) or bending load 
(i = b) , and Q, the shape factor for an ellipse, is given by the square of the 
complete elliptic integral of the second kind (ref. 9) . The half-length of the 
plate h was chosen large enough to have a negligible effect on stress intens- 
ity (h/c = 5) . Values for F, the boundary-correction factor, were calculated 
along the crack front for various combinations of parameters with these ranges; 
0.2 ^ a/c S 2.0, 0.2 ^ a/t ^ 0.8, and 0.2 $ c/b S 0.8. 
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The stress-intensity factors were obtained by use of a nodal-force method, 
the details of which are given in references 6 and 7. In this method, the nodal 
forces normal to the crack plane and ahead of the crack front are used to eval- 
uate the stress-intensity factor. 


RESULTS AND DISCUSSION 

In the following sections, stress-intensity factors for semielliptical sur- 
face cracks in plates subjected to tension or bending loads are presented. The 
effects of crack shape (a/c) , crack size (a/t) , and plate width (c/b) are inves- 
tigated. The stress-intensity factors are compared with other values from the 
literature. Also, a wide-range equation for stress-intensity factors as a func- 
tion of the parametric angle, crack shape, crack size, and plate width is devel- 
oped herein from the present results. 


Tension 

Effects of a/c and a/t .- Figure 4 shows the normalized stress-intensity 
factors as functions of the parametric angle ^ and the ratio of crack depth to 
plate thickness a/t for a semielliptical surface crack (0.2 ^ a/c = 2.0) in a 
plate under uniform- tens ion load. The c/b ratio was taken to be less than or 
equal to 0.2, so that plate width would have a negligible effect on stress 
intensity (less than 1 percent) . For every a/c ratio and parametric angle 
considered, larger values of a/t gave larger normalized stress-intensity fac- 
tors. For low a/c ratios, the maximum normalized stress-intensity factor 
occurs at the point of maximum depth (t = '^/2) , whereas for high a/c ratios 
the maximum factor occurs at the intersection of the crack with the front sur- 
face (4) = 0) . The normalized stress-intensity factors for tension are given in 
table I . 

Effects of c/b .- The normalized stress-intensity factors for a semiellip- 
tical surface crack in a finite-width plate are given in table II for 
a/c = 0.2 or 1.0 and ^ from 0 to tt/ 2 for various values of c/b and a/t. 

In figures 5(a) and 5(b) , the normalized stress-intensity factors for semi- 
circular and semielliptical surface cracks, respectively, are plotted as func- 
tions of 4> for a/t = 0.8 and for c/b from 0.2 to 0.8. For the semicircu- 
lar crack (a/c = 1.0), the stress-intensity factors, which are maximum at 4^ = 0, 
are higher for larger c/b ratios. The dashed lines in figure 5(a) show the 
corresponding stress-intensity factors for a through crack of length 2c. The 
values for a surface crack at 4^ = 0 are about 10 percent lower than the values 
for a through crack of the same length. 

In figure 5(b), the stress-intensity factors for the semielliptical crack 
(a/c = 0.2) are maximum at 4^ = tt/ 2. Again, larger c/b ratios gave higher 
stress-intensity factors. 
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Comparisons with other solutions .- Figures 6(a) and 6(b) show, respectively, 
normalized stress-intensity factors obtained by several investigators for semi- 
circular and semielliptical surface cracks in a finite-thickness plate. The 
results in figure 6(a) are for a/c = 1.0 and a/t = 0.8. The present results 
(solid symbols) are 1 0 to 1 5 percent higher than the results of Smith and Alavi 
(ref. 1) and Kobayashi (ref. 3) for low values of the parametric angle. The 
results of Kathiresan (ref. 4) are fairly close to the present results. All 
solutions agree well near = tt/2. 

Figure 6(b) shows the results for a semielliptical surface crack with 
a/c = 0.2 and a/t = 0.8. The solutions of Smith and Sorensen (ref. 2), 
Kobayashi (ref. 3), and Kathiresan (ref. 4) disagree by 50 to 100 percent. 

Also, these solutions are considerably below the present results. The results 
from reference 2 are generally closer to the present results, although 1 0 to 
25 percent lower. 


Bending 

Effects of a/c and a/t .- Figure 7 shows the normalized stress-intensity 
factors as functions of the parametric angle and a/t ratios for semielliptical 
surface cracks (0.2 ^ a/c = 2.0) subjected to bending loads. (See fig. 3(b).) 
Again, the c/b ratio was less than or equal to 0.2. For all a/c ratios and 
high a/t ratios, the maximum normalized stress-intensity factor occurs at the 
front surface ((}) = 0) . For the lowest a/t ratio considered (0.2), the stress- 
intensity factor variations along the crack front are similar to the tension 
case. Table I gives the normalized stress-intensity factors for these crack 
configurations. 

Effects of c/b .- The normalized stress-intensity factors for a semiellip- 
tical surface crack (a/c = 0.2 or 1.0) are given in table II for various values 
of c/b, a/t, and 4), 

Figures 8(a) and 8(b) show, respectively, the stress-intensity factors for 
semicircular (a/c = 1.0) and semielliptical (a/c = 0.2) surface cracks as func- 
tions of (j) for a/t = 0.8 and c/b frcxn 0.2 to 0.8. For both crack shapes, 
the stress-intensity factor is maximum at the front surface (cj) = 0) . The stress- 
intensity factors are higher for larger c/b ratios. 

Comparisons with other solutions .- Figures 9(a) and 9(b) show stress- 
intensity factors obtained by investigators for semicircular surface cracks in 
a wide plate under bending loads with a/t = 0.6 and 0.8. For a/t = 0.6 and 
4) < tt/ 4 (fig. 9(a)), the present results fall between those of Kobayashi 
(ref. 3) and Kathiresan (ref. 4) . For 4^ > tt/ 4, the results are in close 
agreement. 

The results in figure 9(b) are for a deeper semicircular surface crack 
(a/t = 0.8) than that shown in figure 9(a). Here the present results are 


5 



generally higher by about 10 percent than Kathiresan's and Kobayashi*s results 
for 4) < tt/4. 


Stress-Intensity Factor Equation for the Surface Crack 

For ease of computation, a wide-range equation for the stress-intensity 
factor of a surface crack in a finite plate subjected to tension and bending 
loads has been fitted to the present finite-element results for a/c from 0.2 
to 1.0. To account for the limiting behavior as a/c approaches zero, the 
results of Gross and Srawley (ref. 10) for a single-edge crack have also been 
used. The stress-intensity factor equation for combined tension and bending 
loads is 


Kl 


(St t 


HSb) 



F 


( 2 ) 


for 0 < a/c =1.0, 0 ^ a/t <1.0, c/b < 0.5, and 

approximation for Q, developed by Rawe and used in 


0 ^ (|) S TT. A useful 
reference 11, is 


Q = 


1 + 1 .464( - 

VC 


1 .65 



(3) 


The functions F and H are defined so that the boundary-correction factor 
for tension is equal to F and the boundary-correction factor for bending is 
equal to the product of H and F. The function F was obtained from a sys- 
tematic curve-f itting procedure by using double-series polynomials in terms of 
a/c, a/t, and angular functions of 4). The choice of functions was based on 
engineering judgment. The function F was taken to be 



where 


Ml 


1.13 


0.091 - 

kC 


(5) 
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0.89 

M 2 = -0.54 + ( 6 ) 

a 

0.2 + - 

c 


M 3 = 0.5 


1 .0 

+ 14 

a 

0.65 + - 
c 


.0 



(7) 


g 


+ 


0.1 



sin 4 >) ^ 


( 8 ) 


The function f(j>, an angular function from the embedded elliptical-crack solu- 
tion (ref. 9) , is 


f(J) 



+ sin^ <t> 


V4 


(9) 


The function f^^, a finite-width correction from reference 12, is 


^w 


(TTC a 


sec 


,2bVt 


1/2 


( 10 ) 


The function H, developed herein also by curve fitting and engineering judg- 
ment, has the form 


H = H-| + (H2 - H-j) sinP 


( 11 ) 


where 


p = 0.2 


a 

+ - 


c 


+ 0.6 


a 

t 


( 12 ) 
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For all combinations of parameters investigated and a/t S 0.8, equa- 
tion (2) was within ±5 percent of the finite-element results and the single-edge 
crack solution. (Herein, "percent error" is defined as the difference between 
equation (2) and the finite- element results normalized by the maximum value for 
that particular case. This definition is necessary, especially for the case 
of bending, for which the stress-intensity factor ranges from positive to nega- 
tive along the crack front.) For a/t > 0.8, the accuracy of equation (2) has 
not been established. However, its use in that range appears to be supported 
by estimates based on the concept of an equivalent through crack. Results from 
equation (2) for tension and bending are shown in figures 10 and 11, respec- 
tively, with the stress-intensity factor plotted as a function of 4) for sev- 
eral combinations of a/c and a/t to illustrate the characteristics of the 
equation. 

Equation (2) is used in appendix A to predict the growth patterns of sur- 
face cracks under tension and bending fatigue loads. The predicted growth pat- 
terns are in fair to good agreement with previously published experimental 
measurements made on steel, titanium alloy, and aluminum alloy material. A 
modified form of equation (2) is also used in appendix B to correlate surface- 
crack fracture data from the literature for a brittle epoxy material under 
tension loads. In these data, the ratios of crack depth to plate thickness 
ranged from 0.15 to 1.0 and the ratios of crack depth to crack length ranged 
from 0.3 to 0.84. The equation correlated 95 percent of the data analyzed to 
within ±10 percent of the calculated failure stress. 
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CONCLUDING REMARKS 


A three-dimensional, finite-element stress analysis was used to calculate 
mode I stress-intensity factor variations along the crack front for a wide range 
of semielliptical surface cracks in finite elastic plates subjected to remote 
tension or bending loads. 

For remote tension, the maximum stress-intensity factor occurs at the max- 
imum depth point for small ratios of crack depth to crack length and at the 
intersection of the crack with the front surface for large ratios. For remote 
bending, the maximum value occurs at the intersection of the crack with the 
front surface for ratios of crack depth to plate thickness greater than or equal 
to 0.6, regardless of the ratio of crack depth to crack length. For the lowest 
ratio of crack depth to plate thickness considered (0.2), the stress-intensity 
factor variations along the crack front were similar in the tension and bending 
cases. Also, larger ratios of crack length to plate width give higher stress- 
intensity factors for both tension and bending loads. 

The finite-element results are used to develop a wide-range equation for 
stress-intensity factor for both tension and bending loads. The equation 
applies for any parametric angle, ratios of crack depth to crack length ranging 
from 0 to 1 .0, ratios of crack depth to plate thickness ranging from 0 to 1 .0, 
and ratios of crack length to plate width less than 0.5. For all configurations 
for which ratios of crack depth to plate thickness do not exceed 0.8, the equa- 
tion is within ±5 percent of the finite-element results and the single-edge 
crack solution. For ratios greater than 0.8, no solutions are available for 
direct comparison; however, the equation appears reasonable on the basis of 
engineering estimates. 

The wide-range equation is used in appendix A to predict the growth pat- 
terns of surface cracks under tension and bending fatigue loads. The predicted 
growth patterns were in fair to good agreement with previously published exper- 
imental measurements made on steel, titanium alloy, and aluminum alloy material. 
A modified form of the equation is also used in appendix B to correlate surface- 
crack fracture data from the literature on a brittle epoxy material. In these 
data, the ratios of crack depth to plate thickness ranged from 0.15 to 1.0 and 
the ratios of crack depth to crack length ranged from 0.3 to 0.84. The equa- 
tion correlated 95 percent of the data analyzed to within ±10 percent of the 
calculated failure stress. 

The stress-intensity factor equations presented herein should be useful 
for correlating fatigue-crack-growth rates as well as in computing fracture 
toughness of surface-cracked plates. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
November 19, 1979 
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APPENDIX A 


FATIGUE-CRACK-GROWTH PATTERNS OF SURFACE CRACKS 

The stress-intensity factor equation (eq. (2)) developed for surface cracks 
is used herein to predict fatigue-crack-growth patterns under tension and bend- 
ing fatigue loads. The predicted growth patterns are compared with experimental 
data obtained from the literature for steels, titanium alloys, and aluminum 
alloys . 


Procedure 

The surface-crack configuration considered is shown in figure 1 . Although 
equation ( 2 ) gives the stress-intensity factor at any location along the crack 
front, only the values at the maximum-depth point A and at the front surface B 
were used to predict the crack-growth patterns. (See inset in fig. 1.) The 
cracks were always assumed to be sonielliptical with semiaxes a and c. 

The crack-growth rates were calculated by assuming that the Paris relation- 
ship (ref. 13) between crack-growth rate and stress-intensity factor range is 
obeyed independently at points A and B at the crack front. Thus, 


da 

- = CaAKa'^ 


dN 


(Al) 


dc 

— = CbAKb" (A2) 

dN 


where Ak is the stress-intensity factor range at point A or B, n is an 
exponent to be specified, and and C 3 are the crack-growth coefficients 

for points A and B, respectively. In this paper, n is assumed to be 4, a 
value which has been found to be applicable to a wide range of materials. 
Normally, and Cb are assumed to be equal; however, experimental results 

(refs. 14 and 15) for surface cracks under tension and bending fatigue loads 
show that small semicircular cracks tend to grow semicircular for low a/t 
ratios. Because the stress-intensity factor solution for the small semicircular 
crack (table I) shows that the stress intensity at point B is about 10 percent 
higher than the value at point A, the coefficient €3 was assumed to be 


Cb = c.g^^CA 


(A3) 


SO that a small semicircular crack would be predicted to initially retain its 
shape. Accordingly, equation (A3) was used for all crack configurations con- 
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sidered. On© reason is not equal to Cg may be the changing relationship 

between the stress-intensity factor and the crack-growth rate as the stress 
state changes from plane stress on the front surface to plane strain at the 
maximum-depth point* 

The number of stress cycles required for propagation of a surface crack 
from an initial half-length Cq to a desired half-length Cf was obtained by 
a numerical integration of equation (A2) » This was accomplished by dividing 
the crack extension (Cf - Cq) into a large number of equal increments Ac and 
assuming that each increment was created at a constant crack-growth rate- The 
constant growth rate for each increment was determined from equation (A2) by 
using the crack configuration which existed at the start of that growth incre- 
ment. For each incr^ent of crack advance Ac at the surface, a new increment 
of crack depth Aa was computed from 




^A 

O.9AK3/ 


Ac 


(A4) 


This defined the crack configuration for the next growth increment, and the pro- 
cess was repeated until the crack depth reached the plate thickness. 


Tension 

Figure 12 shows the experimental and predicted fatigue-crack-growth pat- 
terns for surface cracks subjected to tension. The figure shows the a/c ratio 
plotted against the a/t ratio for Ti-6Al-4V titanium alloy (ref. 16), 
9-percent nickel steel (ref. 17), and 221 9-T87 aluminum alloy. 

The experimental procedure for the aluminum and titanium alloys was as 

follows. An electric-discharge machined (EDM) notch was used as a crack 

starter. The vertical bar in figure 12 denotes the range of EDM notch shapes 
(a/c) for the nine aluminum-alloy specimens. Each specimen was subjected to 
constant-amplitude cyclic loading for various numbers of cycles and then stat- 
ically pulled to failure. The data points indicate the final fatigue-crack 
shapes and sizes measured from the broken specimens. Hence, a separate speci- 
men was necessary to obtain each data point. 

The experimental data for the nickel steel were obtained from one specimen 
which was subjected to two- level variable- amplitude loading and then pulled to 
failure. The amplitude change caused "marker bands" to be formed on the crack 

plane. The marker bands in turn were used to define the crack shape and size. 

The dashed curves in figure 12 are the predicted fatigue-crack-growth 
patterns from equation (2) and equations (Al) to (A4) . (Note that n is 
assumed to be equal to 4 and that the growth patterns are independent of the 
magnitude of and Cb* See eq. (A4).) The predicted growth patterns are 

in good agreement with the measurements made on the three materials- The solid 
symbols denote initial crack size and shape. For all initial crack shapes con- 
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sidered, the predicted a/c ratio was about 0.8 when the crack depth became 
equal to the plate thickness. 


Bending 

Figures 13 and 14 show the experimental and predicted fatigue-crack-growth 
patterns for surface cracks in plates subjected to cantilever bending. Although 
the stress-intensity factor equation (eq. (2)) used herein was developed for pure 
bending, the differences between crack-growth patterns for cantilever and pure 
bending are not expected to be large. The a/c ratio is plotted as a function 
of the a/t ratio for aluminum-alloy specimens in figure 13 and steel specimens 
in figure 14. All data points for the 2014-T651 aluminum alloy (ref. 14) were 
obtained from separate specimens, whereas most of the data on the T-1 steel 
(ref. 18) and carbon steel (SS41) (ref. 17) were obtained from single specimens. 
The aluminum-alloy specimens were cycled under constant-amplitude loading and 
then statically pulled to failure, whereas the steel specimens were cycled under 
two- level variable- amplitude loading. The two- level loading produced marker 
bands on the crack plane from which the experimental data were obtained. Again, 
the solid symbols denote the initial crack dimensions. 

In figure 1 3, the predicted fatigue-crack-growth patterns (dashed curves) 
are in good agreement with the experimental data for the aluminum alloy. The 
predicted results show that the cracks tended to approach a common propagation 
pattern, as pointed out by Corn (ref. 14), 

In figure 14, the solid curves are drawn through the T-1 steel data and 
the symbols show the carbon-steel data. Here, the predicted growth patterns 
(dashed curves) are in disagreement with the experimental data for the T-1 
steel. Although the trends are similar, the predicted curves are as much as 
40-percent higher than the experimental curves for T-1 steel, whereas the pre- 
dicted curve for the carbon steel is about 15-percent lower than the experi- 
mental data (symbols) . 

This disagreement between predicted and measured growth patterns for the 
steels prompted a search for additional data. Figure 1 5 shows fatigue-crack- 
growth patterns measured for H-11 steel (ref. 1 5), 4340 steel (ref. 15), D6-AC 
steel (ref. 19), Ti-6A1-4V titanium alloy (ref. 19), 2014-T651 aluminum alloy 
(ref. 14), and the T-1 steel (ref. 18). Again, solid curves were drawn through 
the experimental data. These results show a systematic thickness effect; that 
is, the crack shape (a/c) is generally less for thicker material at a given 
a/t ratio. The upper solid curve (H-11 and 4340 steel) is for the thinnest 
material considered (t = 1 .8 mm), whereas the lowest solid curve (T-1 steel) 
is for the thickest material (t = 25.4 mm). The titanium alloy and the D6-AC 
steel were 6.4 mm thick, the carbon steel was 1 0 mm thick, and the 201 4 -t 6 51 
aluminum alloy was 9.5 mm thick. 

This thickness-related trend may be due to the presence of residual 
stresses in the surface layer (ref. 20) . Tensile residual stresses in the 
surface layer will cause higher stress-intensity factors on the surface than 
factors calculated from equation (2); hence, crack growth will be more rapid 
at the surface. Thus, the a/c ratios will be lower than those predicted for 
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a given value of a/t. In contrast, compressive residual stresses in the sur- 
face layer would cause cracks to grow more slowly on the surface than predicted 
by equation (2) . This would cause higher a/c ratios than those predicted for 
a given value of a/t. 

The dashed curve in figure 15 shows the results predicted by using equa- 
tions (2), (Al), and (A2) . The predicted-results curve was roughly the average 
of all of the experimental data. 
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APPENDIX B 


FRACTURE OF SURFACE- CRACKED BRITTLE MATERIALS 

The stress-intensity factor solutions for surface cracks and for single- 
edge cracks (ref. 10) in plates subjected to tension are used herein to develop 
an equation for fracture of surface-cracked brittle materials. The equation 
is used to correlate fracture data for a wide range of crack shapes and sizes 
in surface-cracked tension specimens made of a brittle epoxy material. 


Stress-Intensity Factor Equation for Fracture 

The application of linear-elastic fracture mechanics to surface-cracked 
specimens is complicated by the fact that the stress-intensity factor solution 
is a function of the parametric angle (t^. For surface cracks with a/c < 0.6, 
the maximum stress-intensity factor occurs at the maximum-depth point, cf) = tt/2; 
for surface cracks with a/c ^ 0.6, the maximum stress-intensity factor occurs 
near the front surface (}) = 0. In fracture analyses, most investigators have 
used the stress-intensity factor at the maximum-depth point (ref. 11). However, 
for a/c ratios greater than about 0.6, some investigators have used the value 
at the front surface because the stress-intensity factor is maximum there. When 
the stress-intensity factor is maximum at the front surface, it may not control 
fracture. The surface material is more nearly in a state of plane stress; hence, 
its resistance to fracture is higher there than in the interior, where nearly 
plane-strain conditions prevail. Thus, for a/c ratios greater than about 0,6, 
fracture may initiate near, though not necessarily at, the front surface. 

In this paper, the maximum stress-intensity factor for a/c < 0.6 and an 
"average” value for 0.6 ^ a/c ^ 1.0 are used to develop an equation for pre- 
dicting fracture. The average stess-intensity factor is the average between the 
values at ({) = 0 and at tt/ 2. For a/c ratios greater than 1.0, an engineer- 
ing estimate similar to that in reference 21 is used. 

For convenience of application, an equation was developed to fit these 
numerical results. The equation is essentially an extension of an equation 
presented in reference 21 . Details of the extension, which involved engineer- 
ing judgment combined with a considerable amount of trial and error, are largely 
omitted. The elastic stress-intensity factor at failure was taken as 


K 


le “ 



(Bl) 


for 0.03 S a/c < °°, 0 ^ a/t < 1.0, and c/b < 0.5. 

Mg, the elastic magnification factor, is 


The remote stress is S^; 
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where p = /tt . The shape factor Q is approximated by 


Q = 1 + 1 .4641 - 


a 

- ^ 1 
c 


Q = 1 + 1 .464U 


The factors M-j and M2 are expressed as 


Ml =1.13- 0.1( - 


0.03 ^ ^ 1 

c 


t; ' * ‘’•“rl 


a 

- > 1 
c 


c/ TT 

M, = , . -U- - , 


f s , 


a 

- > 1 
c 


The finite-width correction f„, obtained from reference 12, is 


= sec 
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For a/c < 0.03, the stress-intensity factor for the single-edge crack plate 
subjected to tension (ref. 10) was assumed to apply, and was written as 


M 


e 



1 . 12 - 



/ a 

21 .71 ^ - 


\t/ 


+ 30.381 - 
t 


(B7) 


Analysis of Fracture Data 

Equation (Bl) has been used to analyze data obtained by Smith (ref. 22) 
in a large number of fracture tests on surface-cracked tension specimens made 
of a brittle epoxy material. (Plane-strain plastic-zone size, based on the 
largest computed fracture toughness, was two orders of magnitude below minimum 
specimen thickness and, hence, consistent with brittle fracture behavior.) 
Thicknesses ranged from 2.5 to 9.5 mm, with 0.15 ^ a/t = 1.0 and 
0.3 ^ a/c ^ 0.84. All specimens were 25 mm wide. 

In this analysis, the specimens were arranged into five groups according 
to their date of manufacture. Specimen thicknesses were constant within each 
group. Fracture of all specimens was assumed to occur at the same value of 
stress-intensity factor (denoted herein by The elastic fracture tough- 
ness for each group of specimens was obtained by averaging the calculated 

stress-intensity factors at failure as 


m 



i=l 


(B8) 


where m is the number of specimens in a group. The values for the five 

groups of specimens were calculated as 677, 682, 713, 723, and 731 kN/m^/^. 

After had been determined, equation (Bl) was used to calculate fail- 
ure stresses. The gross failure stresses were calculated from 


Seal = 



(B9) 


Figure 16 shows the ratio of experimental failure stress S^^p to calculated 
failure stress ^cal Plotted as a function of the a/t ratio. The solid line 
at unity denotes perfect agreement and the dashed lines denote ±1 0-percent scat- 
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ter. The proposed equation cx>rrelated 95 percent of the data analyzed within 
±10 percent for a wide range of a/t and a/c ratios. 
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TABLE I.- NORMALIZED STRESS- INTENSITY FACTORS FOR A SURFACE CRACK IN A 


LARGE PLATE UNDER TENSION OR BENDING LOADS 
[c/b = 0 . 2 ; c/h = 0 . 2 ; V = 0 . 3 ] 



1 .0 

0 

1 .174 


.125 

1 .145 


.250 

1 .105 


.375 

1 .082 


.500 

1 .067 


.625 

1 .058 


.750 

1 .053 


.875 

1 .050 


1 .000 

1 .049 

o 

CM 

0 

0.821 


.125 

.794 


.250 

.740 


,375 

.692 


.500 

.646 


.625 

.599 


.750 

.552 


.875 

.512 


1 .000 

.495 


1 .157 
1 .126 
1.104 
1 .088 
1 .075 
1 .066 
1 ,062 


0.848 
.81 8 
.759 
.708 
.659 
.609 
.560 
.519 
.501 





































TABLE II.- NORMALIZED STRESS-INTENSITY FACTORS FOR A SURFACE CRACK 
IN A FINITE-WIDTH PLATE UNDER TENSION OR BENDING LOADS 


[c/h = 0.2; V = 0.3] 




1 

.153 

1 

.058 

1 

.008 


.968 


.894 


.786 


.660 


.557 


.518 

1 

.574 

1 

.336 

1 

.006 


.721 


.471 


.267 


.113 


.01 5 

- 

.019 

































(a) Specimen model. 


(b) Element pattern on y = 0 



X 

(c) Element pattern around crack front. 


plane. 


Figure 2.- Finite-element model of a plate containing a semielliptical 
surface crack. (r is the radial coordinate.) 
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(a) Semicircular crack. 


Figure 6.- Comparison of calculated stress-intensity factors with results 
from the literature for a deep surface crack in a plate under tension. 
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Figure 12.- Experimental and predicted fatigue-crack growth patterns for a 
surface crack in a plate under tension. (Solid symbols denote initial 
conditions .) 
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Figure 13.- Experimental and predicted fatigue-crack growth patterns for 
a surface crack in an aluminum alloy cantilever plate under bending. 
(Solid symbols denote initial conditions.) 
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